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Abstract 

We derive the effective interaction between two quasiparticles in symmetric nuclear matter resulting 
from the leading-order chiral three-nucleon force. We restrict our study to the L = 0,1 Landau param- 
eters of the central quasiparticle interaction computed to first order. We find that the three-nucleon 
force provides substantial repulsion in the isotropic spin- and isospin-independent component Fq of 
the interaction. This repulsion acts to stabilize nuclear matter against isoscalar density oscillations, a 
feature which is absent in calculations employing low-momentum two-nucleon interactions only. We 
find a rather large uncertainty for the nuclear compression modulus /C due to a sensitive dependence 
on the low-energy constant C3. The eff'ective nucleon mass M* on the Fermi surface, as well as the 
nuclear symmetry energy /3, receive only small corrections from the leading-order chiral three-body 
force. Both the anomalous orbital (/-factor 6gi and the Landau-Migdal parameter (7^jv (characterizing 
the spin-isospin response of nuclear matter) decrease with the addition of three-nucleon correlations. 
In fact, 6gi remains significantly smaller than its value extracted from experimental data, whereas g'j^^j^ 
still compares well with empirical values. The inclusion of the three-nucleon force results in relatively 
small p-wave (L = 1) components of the central quasiparticle interaction, thus suggesting an effective 
interaction of short range. 
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Universe. 



1 



I. INTRODUCTION 

In a recent publication [l| we have studied the quasiparticle interaction in isospin-symmetric 
nuclear matter employing realistic chiral nucleon-nucleon (NN) interactions. The second-order 
calculation performed in ref. [l| explored the order-by-order convergence of the perturbative 
expansion for the quasiparticle interaction as well as the scale dependence resulting from the 
renormalization group evolution of the underlying two-body interaction. Although good 
agreement was found for a number of nuclear observables, such as the effective mass M* at the 
Fermi surface, the nuclear symmetry energy /3 and the spin-isospin response encoded in the 
parameter g'j^j^, both the compression modulus /C of isospin-symmetric nuclear matter and the 
anomalous orbital (jf-factor 6gi were found to differ appreciably from empirical values extrapo- 
lated from collective excitations of finite nuclei. Most seriously, the compression modulus K, of 
nuclear matter (encoded in the Fermi liquid parameter Fq) was found to be negative at both 
first- and second-order across a wide range of cutoff scales. 



Previous studies 



have suggested that interactions induced by the polarization of the 
medium (first studied by Babu and Brown for the case of liquid ^He [91]) provide sufficient 
repulsion to stabilize nuclear matter against isoscalar density fluctuations. In fact, already the 
leading-order contribution to the Babu-Brown induced interaction is quite repulsive jlj], though 
not enough to achieve stability at normal nuclear matter saturation density when employing 



chiral or low-momentum NN interactions. For such potentials, numerous calculations [101413 1 
of the equation of state have revealed the necessity of three-nucleon forces in driving saturation 
toward the empirical density po = 0.16 fm~'^ and energy per particle E = —16 MeV. 

To improve the microscopic description of the quasiparticle interaction, we compute in the 
present work the first-order (perturbative) contribution to the L = 0, 1 Landau parameters 



resulting from the leading-order chiral three-nucleon force. Previous work [14| performed within 
the framework of chiral effective field theory included explicit A-isobar degrees of freedom in 
a calculation of the isotropic central Fermi liquid parameters. In the present study we employ 



the high-precision Idaho N^LO chiral NN interaction 15|] together with the N^LO chiral three- 
nucleon force [l^, which subsumes certain processes involving virtual A-isobar excitations 
considered in ref. [l^. The additional repulsion provided by the leading-order three-nucleon 
force is expected to play an important role for the Landau parameter Fq, but the extent to 
which other nuclear observables respond to additional three-nucleon correlations has not been 
studied systematically in calculations employing high-precision NN potentials. 

The present paper is organized as follows. In Section [TTl we briefly review Landau's theory 



of normal Fermi liquids and discuss the connection between physical observables and various 
Fermi liquid parameters. We then derive analytical expressions for the first-order contribution 
to the L = 0, 1 central Landau parameters arising from the leading-order chiral three-nucleon 
interaction. In section UTTI we present numerical results for the density and scale dependence of 
the quasiparticle interaction supplemented by the second-order contributions of ref. The 
latter were derived from the Idaho N^LO chiral two-nucleon interaction as well as from the 
(universal) low-momentum NN interaction Viow-k- We end with a summary and an outlook. 



II. NUCLEAR QUASIPARTICLE INTERACTION 



A. Landau parameters and nuclear observables 



Fermi liquid theory was introduced by Landau in the 1950's 17| to describe the properties 
of strongly interacting normal many-fermion systems at low temperatures. The low-energy 
excitations about the ground state are long-lived quasiparticles that retain certain features 
of non-interacting (independent) particles but have modified dynamical properties such an 
effective mass and an effective magnetic moment. Fermi liquid theory has been used to describe 
a wide variety of quantum many-body systems, including various types of conductors at low 
temperatures, liquid ^He, nuclear matter, neutron matter and also finite nuclei 
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20|. 



The quasiparticle interaction J^(picriri; P2C2T2) encodes bulk equilibrium and transport prop- 
erties of a Fermi liquid. It arises at second order in the expansion of the energy density in terms 
of powers of the quasiparticle distribution function dn^ar'- 
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where is the single-particle energy and higher-order quasiparticle correlations are neglected. 
The central part of the quasiparticle interaction in spin- and isospin-saturated nuclear matter 
has the form 

J^iPuP2) = f{Pl,P2) + f{Pl,P2)ri ■ T2 + [g{Pl,P2) + 9'{Pl,P2)ri ■ fa] ai ■ (T2 , (2) 

where ai^2 and ri 2 denote the spin and isospin operators of the two nucleons on the Fermi 
surface \pi\ = |p2 | = kp- More generally, the quasiparticle interaction can include (non- 
central) tensor components which modify the stability conditions of nuclear matter 21| and are 
necessary for understanding the magnetic susceptibility 3] and the response of nuclear matter 



to weak probes 



2J]. For two quasiparticles on the Fermi surface, the remaining angular 



dependence of their interaction can be expanded in Legendre polynomials of cos 6 = pi ■ P2'- 

oo 

X{pi,P2) = J2^dkF)PL{cos 9), (3) 

L=0 

where X represents /, /', g, or g', and the angle 6 is related to the relative momentum p = 
^\pi — P2I through the relation p = kp sin {0/2). The coefficients of the expansion in eq.([3]) are 
referred to as the Fermi liquid parameters (FLPs). It is conventional to introduce dimensionless 
FLPs by multiplying by the density of states at the Fermi surface, Nq = 2M*kp /ti"^ where M* 
is the nucleon effective mass and /c^? is the Fermi momentum, leading to 

^ 00 

•^(Pi,P2) = TT E + ^i^i ■ ^2 + {Gl + G^Ti ■ f2)ai • (T2] Pl(cos0). (4) 

For short-range interactions the expansion in (jlj) is typically rapidly converging such that only a 
few constants characterize the dynamics of low-energy excitations. Moreover, individual param- 
eters evaluated at the equilibrium Fermi momentum kp = 1.33 fm~^ are related to properties 
of the quasiparticles and the bulk nuclear medium: 

M* 

Quasiparticle effective mass : — — = 1 + Fi/3, 

Mn 

In-medium orbital gf-factor: oi = ^ H — ^ ] , t^, 

^ ^'2 6(l + Fi/3) 

3fc2 

Compression mo dulus : /C = — ^ (1 + Fo) , 



M* 

Isospin asymmetry energy : (3 = (1 + Fq) 

bpm-isospm response : gj^p^ = — — t^^jq' \^) 



where the anomalous orbital (^-factor Sgi is given by 



and g^^N ~ 13.2 is the strong vrA^ coupling constant. Spin observables, though largely uncon- 
strained experimentally, receive significant contributions from the (non-central) tensor Fermi 



liquid parameters [22 1. 



Recently, we have carried out systematic calculations [l^ of the quasiparticle interaction 
in nuclear matter to second order in many-body perturbation theory employing chiral and 
low-momentum NN interactions. The first- and second-order contributions, shown diagram- 
matically in Fig. [H have the form 

J^''^\piSlti;P2S2t2) = {PlSlti;P2S2t2\V\piSih;p2S2t2) = (12|F|12) (7) 





FIG. 1. First- and second-order diagrammatic contributions to the quasiparticle interaction. The 
wavy hnes symbohze the two-nucleon interaction. 
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where the quantity V = {1 — Pi2)V denotes the antisymmetrized two-body potential and 
rtj = 9{kp — \ kj\) is the usual zero-temperature Fermi distribution. These calculations revealed 
the importance of second-order diagrams in raising the quasiparticle effective mass from M* ~ 
O.TMjv to M* ~ Mjv (both lying in the phenomenological range M*/M = 0.7 - 1.0 [251, bej) as 
well as increasing the symmetry energy from /3 ~ (21 — 25) MeV to /3 ~ (31 — 33) MeV (where 
only the second-order result lies within the range of empirical values (3 = (30—36) MeV 27l.l28[|). 
The isotropic spin-isospin Landau parameter G'q increases at second order, with the effect that 
g'j^j^ changes from a value of around 0.50 — 0.57 to about 0.67 — 0.77, the latter being within the 
range g'j^j^ = 0.6 — 0.7 favored by fits to giant Gamow- Teller resonances in heavy nuclei 29l-l3l[|. 
Despite these encouraging results, the Fermi liquid parameter Fq remained well below —1, giving 
rise to a negative compression modulus K, and a corresponding instability of nuclear matter 
against density fluctuations in the vicinity of saturation density po = 0.16 fm~^. Additionally, 

the anomalous orbital ^f-factor decreased from 6gi ~ 0.3 to 6gi ~ 0.1, which is significantly less 

I, J 

than the empirical value 6gi = 0.20 — 0.26 extracted from giant dipole resonances [32]. This 
feature followed almost entirely from the dramatic increase in the effective mass M* at second 
order. An improved microscopic description of the quasiparticle interaction may require the 
consistent implementation of chiral three-nucleon forces. As a first step in this program, we 
compute here the first-order contribution to the quasiparticle interaction from the N^LO chiral 
three-body force. 

According to eq.([T]) the quasiparticle interaction is obtained by functionally differentiating 



the energy density twice with respect to the quasiparticle distribution functions. For a general 
three-nucleon force, the Hartree-Fock contribution to the energy density is given by 

4n = ^tr.,ntr.,.,tr.3,3 j A^A^i^n,^^ n.^^n,^^ (123|V^37v|123) , (9) 

where V^n = V3Ar(l — P12 — P23 — -P13 + -P12-P23 + -P13-P23) denotes the fully antisymmetrized three- 
nucleon interaction and rig = 6{kF — \kj\) + {27cY 6^ [kj—pj) Suj^jcrjTj- Functionally differentiating 
twice with respect to the two quasiparticle distribution functions then leaves an effective two- 
body interaction containing a single (loop) integral over the filled Fermi sea of nucleons. 

The three-nucleon force employed in the present work is the N^LO chiral three-nucleon inter- 
action 16|, which consists of three components. First, there is a two-pion exchange component 



where qa = 1.29, /^r = 92.4 MeV, m,r = 138 MeV (average pion mass) and qi denotes difference 
between the final and initial momenta of nucleon i. The quantity 

P^l = {-^c,ml + 2c3gl ■ q,) + c^e^^^rldk ■ (gl x g,) (11) 

involves three terms proportional to the low-energy constants Ci,C3 and C4, respectively. The 
summation 'Yl,i^j^k runs over the six permutations of three nucleons. The one-pion exchange 
component of the three-nucleon interaction is proportional to the low-energy constant Cd'- 

[1.) _ sr^ QACp aj-qj 

and finally the three-nucleon contact interaction introduces the low-energy constant ce'- 

^^^'= E^^-n-' (13) 

where = 700 MeV sets a natural scale. 

The low-energy constants of V^f^^ appear already in the two-pion exchange contribution to 
the nucleon-nucleon interaction and can therefore be fit to NN scattering phase shifts. The anal- 



ysis of the Nijmegen group [33| resulted in the values Ci = —0.76 GeV S C3 = —4.78 GeV \ and 



C4 = 3.96 GeV"\ while Entem and Machleidt [l5| obtain ci = -0.81 GeV-\ C3 = -3.2 GeV~\ 
and C4 = 5.4 GeV^^. The low-energy constants cd and ce of the mid- and short-range chiral 
three-nucleon interaction can be constrained by properties of nuclear few-body systems. In the 
present work we employ two versions oi cr, and ce obtained by fitting the binding energies of 
A = 3 nuclei together with the half-life of with the result ^341] : 

Co = -0.20, ce = -0.205, (14) 
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FIG. 2. Diagrammatic contributions to the quasiparticle interaction in symmetric nuclear matter gen- 
erated from the two-pion exchange three-nucleon force. The short double-Une symboUzes summation 
over the filled Fermi sea of nucleons. Reflected diagrams of (2) and (3) are not shown. 



or by fitting the binding energies of '^H and "^He Uj : 

CD = -2.062, ce = -0.625. (15) 

In the first set the leading-order three-body force (with 01,3,4 coefficients of Entem and Mach- 
leidt) was used together with the Idaho N^LO chiral two-nucleon interaction, while in the 
second set the authors employed the Nijmegen low-energy constants 01,3,4 together with the 
low-momentum NN interaction Mow-k at the resolution scale A = 2.1 fm~^. Employing these 
two versions of the chiral three-nucleon force (combined with two-nucleon interactions at dif- 
ferent scales) provides a means for assessing theoretical errors at this order in the perturbative 
expansion. 



B. Diagrammatic calculation 

We begin by considering the isotropic {L = 0) Fermi liquid parameters of the central quasi- 
jarticle interaction. They are obtained by angle- averaging the (in-medium) effective interaction 

Mkf) = j dn.dn, {piJ2\v^n\pi.P2) 

= foikf) + go{kf) CTi • as + foikf) n ■ T2 + g'^ikf) Bi • 0^2 n ■ , (16) 

where IP1P2) denotes an antisymmetrized two-nucleon state. Both quasiparticle momenta pi 
and P2 he on the Fermi surface so that |pi| = |p2| = kj. There are six topologically distinct 
one-loop diagrams, shown in Figs. |2] and [3l contributing to the effective interaction VJ^^. Using 
the abbreviations cr = Bi ■ a2 and r = ri ■ 72, we find for the crossed term of diagram (1) in 
Fig. [2] the following contribution 

^.(fc,)'--'" = (3 - <.)(3 - r)0jj\ ^'''l-^f - c,u' + (c, - cO| ln(l + 4„=)} , (17) 
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FIG. 3. Diagrammatic contributions to the quasiparticle interaction in symmetric nuclear matter 
generated by one-pion exchange and contact three-nucleon forces. Reflected diagrams of (4) and (5) 
are not shown. The last two diagrams (5) and (6) contribute only to the L = Landau parameters. 

with u = kp/mT^. In this diagram, the two pions carry equal momenta and therefore the 
term in eq.ffTT]) proportional to C4 does not contribute. Moreover, the direct term (with zero 
momentum transfer) of this pion self-energy correction vanishes trivially. Analogously, only the 
crossed term from the one-pion exchange vertex correction (diagram (2) in Fig. [2]) is nonzero 
with the analytical result: 



Mkf) 



(med,2) 



(3-^)(3-t)- 



g%m^^ j 3ci 



X 



(247r)2/n8w^ 
+ Au^ - (1 + 4^2) ln(l + Au^ 

2 



- ln(l + 4u^ 



+ 



+ C4 



— (^4m^ - ln(l + 4m^) j arctan 2u + 
ln(l + 4^2) + 



+ C3 

12u* - 16u^ - 30^2 - 9 



64^? + 



20^3 1 3 



24m5 



ln(l + 4u^) - 4m^ ) arctan 2u 



+ - + ^ „ 



3 + 16^2 - 48u^ 
64^ 



ln^(l + 4m^ 



80m6 + 12m^ - 30^2 - 9 , ^ 2^ 



28u^ 1 3 ■ 

—- + 13U + - + — 
3 u Au'^ 



The Pauli-blocked two-pion exchange component (diagram (3) in Fig. [2]) has both a nonvan- 
ishing direct and crossed term. Their sum takes the form 



(med,3) 



3 

24(c3 - ci)u - 8c3M^ + (3c3 - 4ci)- ln(l + 4^^^ 

u 



+ 6(6ci - 5C3) arctan 2u + (3 - o-)(3 - r) ^ /" (Y^ - X^) 
+ (1 + o-)(1 + t) 



dx 



3c, z2 + ^(X2 + 2F2) 



(19) 



where the auxihary functions X, Y, and Z arising from Fermi sphere integrals over a pion 
propagator read 

, 5x2-3^2-3 Ax^ + 3{1 + - x^)\ 1 + {u + 
= A-x + W 1 + iu-xr ' ^''^ 

ZKx) = l + ^^^^ln ; + i- + t . (22) 

^ ^ 4mx 1 + (u - x)2 ^ ^ 

Note that the direct term contributes only to the spin- and isospin-independent Landau pa- 
rameter Fq. 

There are two diagrammatic contributions from the mid-range one-pion exchange chiral 
three-nucleon force, labeled as (4) and (5) in Fig. 121 The crossed term from the Ivr-exchange 
vertex correction (diagram (4)) leads to the contribution 

MkfY-'"''^ = (3 - ^)(3 - r) //^7/ Uu' - nln(l + Au^)] , (23) 

(247r)^/4A^ L J 

and the sum of direct and crossed terms from diagram (5) yields 
:Fo{kfY-'^'^ = {3-a-r " (|^^{^ - ^ + -^^^^2^ - i- ln(l + 4^^)} . (24) 

Finally, the three-nucleon contact term generates a contribution proportional to the nuclear 
density p = 2kp/37i'^. The sum of direct and crossed terms reads 

Mkfr-'''' = {cr + r + ar- 3)^^ • (25) 

We note that for all three-body contributions the spin-spin and isospin-isospin components of 
the quasiparticle interaction are equal, go = /q. 

The p-wave (L = 1) Fermi liquid parameters are obtained by weighting the angular integrals 
by the first Legendre polynomial 3Pi(cos^^) = 3pi ■ p2'- 

J'likf) = ^/ d^^dn, (pi ■P2){puP2\Vnn\PuP2) 

= fi{kf) + 9i{kf) Si-a2 + f'i{kf) Ti • fs + g'i{kf) ai • (T2 n ■ . (26) 

After this weighting, the two short-range contributions (diagrams (5) and (6) in Fig. [3]) to 
the quasiparticle interaction vanish. We provide the expressions for the four remaining pieces 
below. The crossed term from the pion self-energy correction takes the form 

+ [2(C3 - ci)'t. + + in')\ . (27) 



The crossed term from the one-pion exchange vertex correction reads 



(med,2) 
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X 
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+ 



(167r)2/n4^i' 
Su"^ + Au^ - (1 + Au^) ln(l + Au^) 

_C3_ 



- (1 + 2^2) ln(l + 



4 (^4^2 - (1 + 2m^) ln(l + 4^2)) arctan 2m 



r ln^(l + Au^) 



56^6 - 32m^ - 60m4 - 48^2 - 9 , , , 

+ ^ + ' 

+ 1!^(7-4m2)-14m'' + 6m+ — 

3 2u 
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+ 



8 (^(1 + 2^2) ln(l + Au^) - 4^2^ arctan 2m 



C4 

3m^ 

3 + 22^2 - 16m^ - 96^6 



ln2(l + 4^2 



32^5 

160m^ - 136^6 - 60m^ - 48^2 - 9 , , on 

+ + ' 

- —{4u^ + 17) + 34^3 + 5u+^ 
3 2it 

The crossed term from Pauh-blocked two-pion exchange is given by 



ct)(3-t)^(3X,2 + 2X,2-3X2) 



+ (1 + o-)(l + r) [3ci(Z„^ + 2^2) + 1(3X2 ^ g^2 ^ 



with auxihary functions 



^ , , a; ti^ - - 1 1 + + x)2 
Za{u, x) = - H In — — , 

Zi,{u, x) = : + ^ In ■ ^ ^ 



4^2; 



1 + (u- xf 



X ( )-- - + + 1 In 1 + (m + 
" ' u Av?x 1 + (rt — x)^ ' 



1 



2x^(1 +M2)-3a;^- 3(1 + M 



2\2 



+ 



1 + ?7,2 + X^ 



32u^x^ 
1 



3(l + M^-x^)^ + 4x 



In 



Xc{u, x) = - — ^ 3(1 + ii^)2 - %u^x^ - 2>x 



+ 



8^0:^ 
3 

32^2^3 



1 + (m + x)2 

1 + (m - X)2 ' 

l + (M + a;)2 
1 + (m - x)2 



(x^ - 1 - M^) (1 + + a;2)2 - 4mV 



In 



(28) 



(29) 



(30) 
(31) 
(32) 

(33) 

(34) 



Finally, the only nonvanishing term from the mid-range three-nucleon force is the crossed term 
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from diagram (4) in Fig. |3l which reads 



^,(,^,,....., . (3 _ .)(3 _ r)5|g^{4,. - (2. + i) ,„(1 + 4„=)} . 



(35) 



Again, one observes that for all contributions, f[ = gi- 

This completes the first-order calculation of the L = 0, 1 Landau parameters arising from 
the N^LO chiral three-nucleon force. A good check of the formulas in eqs.(17)-(35) is given by 
their Taylor expansions in fci;'. The leading terms are of the form: CEkp,CDkp,cikp,C3^4^kp. Note 
that full consistency with the second-order calculation in ref. |l| would require the inclusion of 
the subleading N^LO chiral three-body force (recently derived in ref. j36|). Investigations along 
these lines are in progress. 

III. RESULTS 

In this section we study the density-dependence of the L = 0, 1 Landau parameters derived 
in the previous section from the leading-order chiral three-nucleon interaction. The values 
of the five low-energy constants occuring at this order have significant uncertainties, and we 
employ two different sets which have been fit to reproduceproperties of A = 3, 4 nuclei. These 
contributions are then combined with the results of ref. \M\ for the second-order quasiparticle 
interaction computed with the corresponding two-nucleon interactions. The impact on various 
nuclear observables is discussed. 

In Fig. IHwe plot the L = 0, 1 Fermi liquid parameters (in units of fm^) as a function of the 
nuclear density (normalized to that of saturated nuclear matter po = 0.16 fm~^) employing the 
set of low-energy constants 

ci = -0.81GeV"\ C3 = -3.2GeV-\ C4 = 5.4GeV"\ cd = -0.20, ce = -0.205, (36) 

which have been used with the Idaho N^LO chiral two-nucleon interaction to reproduce the 
binding energies of A = 3 nuclei and the triton half-life. We observe that for densities greater 
than p ~ Po/4, both the L = and L = 1 Fermi liquid parameters depend approximately 
linearly on the density, which is not immediately obvious from the analytical expressions given 
in the previous section. The largest effect on the quasiparticle interaction is a dramatic in- 
crease in the isotropic spin- and isospin-independent component, /o, which at nuclear matter 
saturation density is about four times larger than any of the other contributions to the L = 
parameters. In fact, both qq and /g remain negligibly small for all densities considered here, 
while the parameter q'q decreases monotonically with the density and reaches the value — 0.3fm^ 
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FIG. 4. Density-dependent L = 0,1 Fermi liquid parameters obtained from the leading-order chiral 
three- nucleon force with low-energy constants given in eq. (I36p . 



at nuclear matter saturation density po- The L = 1 Fermi liquid parameters all decrease with 
increasing density. The relatively small change in /i provides only a small downward correction 
to the quasiparticle effective mass M*. This observation, combined with the small change in 
/g discussed above, indicate that three nucleon forces have a relatively small effect on the nu- 
clear symmetry energy /3. The significant reduction in the parameter /( decreases the nucleon 
anomalous orbital gf-factor 6gi considerably below its empirical value. 

In Table [I we show the various contributions to the L = 0, 1 Fermi liquid parameters arising 
from the six topologically-distinct diagrams of Figs. |2] and [31 These are labeled as V^^' * for 
i = 1, ... ,6. Contributions from the long-range two-pion exchange component of the three- 
body force are significantly larger than those from the mid- and short-range three-body forces. 
However, since both V^^'^ and V^^''^ are medium modifications to one-pion exchange but 
enter with opposite sign, to a large extent they cancel in all Fermi liquid parameters. This leaves 
the Pauli-blocked two-pion exchange diagram V^^^*^'^ as the dominant contribution, except in 
the Fermi liquid parameters go = /q and /i where its effects are surprisingly small. It is 
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kp = 1.33 fm-i 




Turned, 1 


Turned, 2 


T rined, 3 


Turned, 4 


T rined, 5 


Turned, 6 


Total 


/o [fm2] 


0.976 


-1.196 


1.277 


-0.018 


-0.036 


0.215 


1.218 




-0.325 


0.399 


-0.011 


0.006 


0.012 


-0.072 


0.009 




0.108 


-0.133 


-0.209 


-0.002 


0.012 


-0.072 


-0.295 


/i [fm2] 


-0.359 


0.216 


0.064 


0.005 








-0.073 


51 = /( [fm^] 


0.120 


-0.072 


-0.278 


-0.002 








-0.232 




-0.040 


0.024 


-0.164 


0.001 








-0.179 



TABLE I. Individual contributions to the L = 0, 1 Fermi liquid parameters arising from the six 
diagrams in Figs. [2] and [3l The contributions are evaluated at nuclear matter saturation density 
po = 0.16 fm^^ corresponding to a Fermi momentum of kp = 1.33 fm^^. The values of the low-energy 
constants of the leading-order chiral three- nucleon force are those given in eq. (j36p . 



instructive to compare the final results to what would be obtained in the chiral limit (m^r — > 0): 

{hY = 1.287fm2, (^o)^ = {fX = 0.136 fm^, {g',Y = -0.424fm2 
[hY = -0.073 fm^, {g^Y = UiY = -0.268 fm^, {g[Y = -0.203 fm^, (37) 



which are not very different from those at the physical pion mass. 

To assess the scale dependence of the Fermi liquid parameters resulting from the choice of 
momentum cutoff in the nuclear interaction, we consider next the set of low-energy constants 
employed in ref. ll| to fit A = 3,4 binding energies using in addition a low-momentum NN 
interaction Viow-k defined at the resolution scale A = 2.1 fm~^: 

ci = -0.76GeV"\ cg = -4.78 GeV~\ C4 = 3.96GeV-\ c^^ = -2.06, c^j = -0.625. (38) 

In Fig. [5] we plot the density-dependent Fermi liquid parameters resulting from the first-order 
perturbative contribution from the chiral three-nucleon force. Qualitatively, many of the trends 
observed with the previous set of low-energy constants (eq. (136|) ) remain. However, the value 
of the isotropic spin-independent and isospin-independent Fermi liquid parameter /o is signifi- 
cantly larger. At nuclear matter saturation density, its value of /o = 2.49 fm^ is twice as large 
as that observed in Fig. |H This repulsion more than compensates for the additional attrac- 
tion that results from employing low- momentum two-nucleon interactions [l|. The parameters 
go and /g, which were nearly zero across the wide range of densities shown in Fig. HJ now 
achieve a value of — 0.24fm^ at saturation density. All of the L = 1 Fermi liquid parameters 
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FIG. 5. Density-dependent L = 0,1 Fermi liquid parameters obtained from the leading-order chiral 
three- nucleon force with low-energy constants given in eq. (I38p . 



exhibit a similar density dependence, and at saturation density lie within the range — 0.2fm^ to 
— 0.3fm^. The Fermi liquid parameter /i, determining the quasiparticle effective mass M*, is 
twice as large as that obtained using the low-energy constants given in eq. (!36l) and thus gives 
rise to a more visible reduction of M*. 

We display in Table Hi] the decomposition of the Fermi liquid parameters into contributions 
from the six unique topologies show in Figs. [2] and [31 Although the mid- and short-range 
components of the three-nucleon force are now larger in magnitude, taken together they largely 
cancel (in the L = FLPs) and again give a relatively small contribution to the total sum. 
The largest differences between the values given in Tables HI and HTl arise from the parameter C3. 
Employing the low-energy constants in eq. (!36ll except with C3 = —4.78 would yield the Fermi 
liquid parameters 

/o = 2.237fm2, (?o = /o = -0.150 fm^, 
/i = -0.210 fm^ ^1 = /{ = -0.293 fm^ 

which are very close to those in Table [Tll 

Finally, we combine these results with previous calculations [l] in which the second-order 
contribution to the quasiparticle interaction was computed to second order employing chiral 



9o 



-0.341 fm^ 
-0.265 fm^ 



(39) 



14 



kp = 1.33 fm-i 




Turned, 1 


Turned, 2 


T rined, 3 


Turned, 4 


T rined, 5 


Turned, 6 


Total 


/o [fm2] 


1.412 


-1.171 


2.149 


-0.189 


-0.369 


0.655 


2.488 




-0.471 


0.390 


-0.124 


0.063 


0.123 


-0.218 


-0.237 




0.157 


-0.130 


-0.270 


-0.021 


0.123 


-0.218 


-0.359 


/i [fm2] 


-0.570 


0.394 


-0.082 


0.054 








-0.203 


51 = /( [fm^] 


0.190 


-0.131 


-0.333 


-0.018 








-0.293 




-0.063 


0.044 


-0.250 


0.006 








-0.263 



TABLE II. Individual contributions to the L = 0, 1 Fermi liquid parameters arising from the six 
diagrams in Figs. [2] and [3l The contributions are evaluated at nuclear matter saturation density 
po = 0.16 fm^^ corresponding to a Fermi momentum of kp = 1.33 fm^^. The values of the low-energy 
constants of the leading-order chiral three- nucleon force are those given in eq. (j38p . 



and low-momentum NN interactions. The density dependence of the dimensionful Fermi liquid 
parameters obtained from the low-momentum NN interaction Viow-k at a cutoff scale A = 
2.1 fm~^ together with the associated three-nucleon force is shown in Fig. |6] (qualitatively similar 
results hold also for the unevolved NN interaction considered in this work). Both the L = 0, 1 
components of the spin- and isospin-independent quasiparticle interaction vary strongly with the 
density, and the inclusion of three-nucleon forces significantly enhances the trend observed in the 



second-order calculation including two-body 



'orces only [l| (note that in contrast to the current 



calculation, the results shown in Fig. 3 of ref. l| did not yet include self-consistent single-particle 
energies in the denominators of the second-order terms). The parameters associated with the 
spin- and isospin-dependent components of the quasiparticle interaction, on the other hand, 
are much more stable to variations in the density (in qualitative agreement with observations 
made in ref. jl^). In fact, for these components the three L = terms have a remarkably 
similar density dependence, and the three L = 1 terms exhibit this same property. 

We now discuss the implications for physical observables. We display in Table IIIII the 
dimensionless Fermi liquid parameters (obtained by multiplying the dimensionful parameters 
by the density of states at the Fermi surface A^o = 2M*kF/'n''^) and the related quasiparticle and 
bulk nuclear matter properties resulting from calculations including two-nucleon interactions 
alone as well as including the leading-order chiral three-nucleon force. The former are denoted 
with superscript 2N while the latter are denoted with 3A^. The most important feature is 
a qualitative change of the isotropic spin- and isospin-independent Fermi liquid parameter 
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Fq, which was large and negative without three-body forces but which now attains a mean 
value that is positive and relatively small. In fact, the compression modulus /C computed 



with the low-momentum interaction is now quite large (/C ~ 530 MeV) compared to va' 
}C = (200 — 300) MeV extracted from studies of giant monopole resonances in heavy nuclei 



ues 



37- 



. Given that the use of low-momentum two- and three-body interactions yields a quite good 



description of nuclear matter saturation already in (approximate) second order perturbative 



calculations 



11 



13| . it is likely that the large value of Fq obtained with \4ow-k in the present 
work results from the absence of second-order contributions from three-body forces. Aside 
from the large uncertainty in the value of Fq, nearly all the other Fermi liquid parameters now 
exhibit much less scale dependence. This feature has also been observed in previous calculations 
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, isl , where it was found that including three-nucleon forces and higher-order corrections in 
perturbation theory dramatically reduced the scale dependence of the nuclear matter equation 
of state. The effective mass M* and the nuclear symmetry energy /3 receive relatively small 
corrections from three-body forces at this order and to within errors remain in agreement 
with their empirical values. The observable that is not well described in this calculation is 
the anomalous orbital ^f-factor 6gi which remains considerably below the value extracted from 
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Go 


F' 




Fi 


Gi 


F[ 




M*/Mn 


JC [MeV] 


(3 [MeV] 




y2N 


-1.64 


0.35 


1.39 


1.59 


-0.13 


0.50 


0.58 


0.47 


0.96 


-150 


31 


0.12 


y2N 

low— k 


-1.98 


0.58 


1.94 


2.14 


0.38 


0.83 


0.87 


0.80 


1.13 


-190 


32 


0.07 


T/3iV 


-0.15 


0.35 


1.36 


1.19 


-0.22 


0.21 


0.28 


0.24 


0.93 


200 


31 


0.09 


T/3JV 
low— k 


1.48 


0.22 


1.45 


1.48 


0.08 


0.37 


0.41 


0.39 


1.03 


530 


29 


0.05 



TABLE III. Dimensionless Fermi liquid parameters obtained by summing the second-order contribu- 
tion from the Idaho N'^LO potential and low-momentum NN interactions V|ow-k(A = 2.1 fm~^) at 
kp = 1.33 fm~^. Hartree-Fock self-energy insertions, as parameterized in ref. [It], are included in the 
second-order diagrams. 



giant dipole resonances of nuclei. This quantity, however, is particularly sensitive to the value 
of Fi and the quasiparticle effective mass (see eq.Q), and therefore even a moderate decrease 
of the nucleon effective mass M* due to higher-order perturbative corrections could remedy 
the situation. We summarize below the results for selected nuclear observables, including 
also theoretical uncertainties estimated by choosing different resolution scales and low-energy 
constants: 



Effective mass : 



M* 



M, 



N 



0.98 ±0.05, 
0.07 ±0.02, 



Anomalous orbital (yf-factor: Sgi = 

Compression modulus : K, = (370 ± 160) MeV, 
Isospin asymmetry energy : /3 = (30 ± 1) MeV, 
Spin-isospin response : g'j^^ = 0.55 ± 0.03 



[0.7-1.0] 
[0.20 - 0.26] 
[200 - 300] MeV 
[30 - 36] MeV 
[0.6-0.7], 



(40) 



where the quantities in the last column are estimates of the empirical values (see Section UTIA). 
A final interesting feature of the three-nucleon force is that the contributions to the L = 1 
Fermi liquid parameters are all negative, with the result that the quasiparticle interaction at 
the Fermi surface appears to be short ranged. 



IV. CONCLUSION 

In this work we have performed the first systematic inclusion of the N^LO chiral three- 
nucleon force in calculations of the quasiparticle interaction in nuclear matter. Analytical 
expressions for the L = 0, 1 Fermi liquid parameters were derived, and numerical results have 
been extracted for two different sets of low-energy constants (c^, cd and ci^s^) corresponding to 
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two choices of resolution scale A = (500,414) MeV for the two- and three-nucleon interactions. 
In comparison to calculations employing two-nucleon forces only, we find a significantly smaller 
dependence of the Fermi liquid parameters on the choice of this cutoff scale. At this order in 
the perturbative calculation, the bulk properties of nuclear matter and the dynamical features 
of quasiparticles are in general well described. In particular, the compression modulus /C of 
nuclear matter attains a reasonable value only with the inclusion of three-nucleon correlations, 
which introduce large additional repulsion in the effective interaction. Although an analogous 
effect has already been observed in calculations of the nuclear matter equation of state, the fact 
that the three-body force gives a disproportionately large contribution to Fq, compared to all 
other Fermi liquid parameters, comes as a surprise. The remaining contributions to the L = 0, 1 
Fermi liquid parameters from the chiral three-nucleon force are (essentially) all negative. This 
results in diminished L = 1 Landau parameters and an effective interaction of apparent short 
range. A remaining unresolved problem concerns the description of the anomalous orbital g- 
factor 6gi (which arises from interactions with the nuclear medium). In the present work we 
find corrections to free space orbital (yf-factor that are much too small, a feature which mainly 
results from a large nucleon effective mass (M* ~ Mn). The role played by three-nucleon 
force contributions beyond the Hartree-Fock approximation and the effect of the subleading 
(N^LO) chiral three-body force 36| remain topics of interest for future calculations. A possible 
strategy is to represent three-nucleon force contributions in terms of density-dependent two- 
body forces 35| and to employ these in second-order calculations. We have verified that treating 
three-nucleon force contributions in this approximation gives good agreement with their exact 
evaluation at first order. 



We thank A. Schwenk for helpful discussions. 
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